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ABSTRACT 

We present 20 band photometry from the far-IR to X-ray in the Spitzer IRAC dark field. The bias 
for the near-IR camera on Spitzer is calibrated by observing a ~ 20' diameter "dark" field near the 
north ecliptic pole roughly every two-to-three weeks throughout the mission duration of Spitzer. The 
field is unique for its extreme depth, low background, high quality imaging, time-series information, 
and accompanying photometry including data taken with Akari, Palomar, MMT, KPNO, Hubble, 
and Chandra. This serendipitous survey contains the deepest mid-IR data taken to date. This 
dataset is well suited for studies of intermediate redshift galaxy clusters, high redshift galaxies, the 
first generation of stars, and the lowest mass brown dwarfs, among others. This paper provides a 
summary of the data characteristics and catalog generation from all bands collected to date as well as 
a discussion of photometric redshifts and initial and expected science results and goals. To illustrate 
the scientific potential of this unique dataset, we also present here IRAC color color diagrams. 
Subject headings: galaxies: photometry — cosmology: observations 



1. INTRODUCTION 

Deep, multi-wavelength surveys are an important tool 
in studying the formation and evolution of galaxies. In- 
frared data is particularly useful in these surveys. Near- 
IR wavelengths provide a direct measurement of the stel- 
lar content of galaxies in absence of a dominant AGN 
contribution at redshifts as high as three. The longer 
wavelengths sample emission primarily from polycyclic 
aromatic hydrocarbons (PAHs), as well as direct ther- 
mal emission from hot dust. These are markers of star 
forming galaxies and active galactic nuclei (AGN). 

The dark field is a multi-wavelength survey based 
on observations carried out as part of t he calibra- 
tion for the Infrared Array Camera (IRA C: | Fazio et al 
2004) on the Spitzer Space Telescope (jWerner et al 
2004). Rather than being designed around specific 
science objectives, like many other Great Observatory 
surveys (H DF, CDF, SWIRE, GOODS, GEMS COS- 
MOS, etc.; IWilliams et all H9M IGiacconi et all |200H 
Lonsdale et al.1 120031: iGiavalisco et all 120041 : iRix et all 
20041: IScoville et all I2007n ~ the survey area and obser- 



vation strategy was selected based on the calibration re- 
quirements for IRAC. Observations were taken to mea- 
sure the dark current and biases of the four detector 
arrays (see £12.1. ip near the beginning and end of each 
IRAC observing period (or "campaign"), resulting in 
2-5 dark field measurements every two-to-three weeks 
throughout the entire duration of the mission. Each of 
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these observations in the dark field collects roughly two 
hours of integration time total in each of the four bands in 
a single 5' x 5' field of view. Each new visit is observed in 
a slightly different position due to spacecraft orientation 
and rotation. The center of the field is located at 17 h 40 m , 
+69 d (J2000), very near the north ecliptic pole (NEP). 
Over the course of the mission, the observations have 
filled in a region roughly 20' in diameter. This has cre- 
ated the deepest ever mid-IR survey, exceeding the depth 
of the deepest planned regular Spitzer surveys over sev- 
eral times their area. Furthermore, this is the only field 
for which a 5+ year baseline of m id-IR periodic observa- 
tions exists(see iFrost et al.ll200"9l for more inforrmation 
on the time-series data). 

The resulting observations are unique in several ways. 
The dark field lies in the lowest possible region of zodia- 
cal background, the primary contributor to the infrared 
background at these wavelengths, and as such is in the re- 
gion where the greatest sensitivity can be achieved in the 
least amount of time. The area was also chosen specifi- 
cally to be free of bright stars and very extended galax- 
ies, which allows clean imaging to very great depth. The 
field lies in the northern continuous viewing zone, a 10 
degree radius region centered on the NEP that is always 
visible to Spitzer, which means many opportunities ex- 
ist to re-observe the field essentially continuously. The 
observations are done at many position angles (which 
are a function of time of observation) leading to a more 
uniform, circularized, final point spread function (PSF). 
Finally, because the calibration data are taken directly 
after anneals, they are more free of electronic artifacts 
such as hot pixels and latent images than ordinary guest 
observer (GO) observations. 

To enhance the IRAC dark field dataset we have ob- 
tained imaging data in 16 other bands including Palo- 
mar/MMT/KPNO u', g', r', F, z', J, H, K, HST /ACS 
F814W, Akari 4, 11, 15, &18^m, Spitzer MIPS 24&70/im, 
and Chandra ACIS-I. Although the entire dark field is 
> 20' in diameter, because of spacecraft dynamics the 
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central ~ 15' is significantly deeper and freer of artifacts. 
Therefore, it is this area which we have matched with 
the additional observations. All space-based datasets are 
publicly available through their respective archives. 

This paper is structured in the following manner. In 
fj2] we present the data for all 20 bands, its reduction 
as well as source extraction and photometry. Details of 
the composite catalog are presented in Sj3] Photometric 
rcdshifts are derived and errors calculated in J|J Lastly, 
a small discussion of the science goals are presented in 
S}5] Throughout this paper we use Ho — 70km/s/Mpc, 
= 0.3, Qa = 0.7. All photometry is quoted in the 
AB magnitude system. 

2. DATA 

The details of our dataset are summarized in Table [T] 
and in the sections below. Data are presented roughly in 
order of decreasing wavelength with the exception that 
the IRAC data are discussed first. Figure Q] shows the 
filter curves of the optical through infrared bands dis- 
cussed below. Figure [2] shows a representative sample of 
data thumbnails. 

2.1. Spitzer Space Telescope 
2.1.1. IRAC 

IRAC simultaneously observes with four detectors at 
3/6. 4/5. 5/8. & 8/0 microns. Each detector has a 5.2' x 
5.2' field of view with 1'.'2 pixels. Images are acquired in 
pairs (3.6 and 5.8/im, 4.5 and 8.0/zm) of two adjacent 
fields, separated by 6.5', using beam splitters. 

This work is based on a preliminary combination of 
calibration data taken from Dec 2003 through Jan 2005; 
«30% of the expected dark field depth not including the 
warm mission, see Figure [3] The 3.6, 4.5, 5.8, and 8.0 
/im data have a total exposure time of 40 hours per filter 
covering 380 square arcminutes to a limiting depth of 0.2, 
0.17, 0.11,0.11 fi3y respectively. The maximum per pixel 
exposure time is 28 hours at 3.6 and 5.8 /im and 21 hours 
at 4.5 and 8.0/^m. The median exposure time per pixel is 
8-10 hours; see below. Raw astrometry from the Spitzer 
pointing control system, based on inputs from the star- 
tracker and gyroscopes, is used as the input to a pointing 
refinement program. The absolute positions of stars from 
the Two Micron All Sky Survey (2MASS) along with the 
relative positions of stars common to different individual 
IRAC frames are then used to improve the astrometry 
from the raw spacecraft astrometry. Final astrometric 
accuracy for an individual frame is ~ 0.3" rms. Because 
field rotation during the mission should average out many 
systematic effects, the final astrometric accuracy of the 
IRAC dark field mosaic is expected to be better than 
this. The Basic Calibrated Data (BCD) products pro- 
duced by the Spitzer Science Center (pipeline versions 
Sll, S12, S13) were further reduced using a modified ver- 
sion of the pipeline developed for the Spitz er Wide area 
Infra red Extragalactic survey (SWIRE) ((S urace et aLl 
2005). This pipeline primarily corrects image artifacts 
and forces the images onto a constant background (ne- 
cessitated by the continuously changing zodiacal back- 
ground as seen from Spitzer). The data were coadded 
onto a regularized 0.6" grid using the MOPEX software 
developed by the Spitzer Science Ce nter (SSC) as de- 
scribed in Mak ovoz fc Marleaul (|2005l ). Experiments we 
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Fig. 1. — Filter transmission curves for the optical to far-IR 
bands. All curves are normalized for comparison. The upper plot 
shows optical through near-IR. Of the near-IR bands, only the 
Palomar WIRC J, H, and K curves are shown for clarity. The 
lower plot shows the Spitzer and Akari (dashed lines) mid-IR to 
far-IR bands. 




Fig. 2. — Representative images of the field, 3 arcminutes on a 
side. The order of the images is u' , r' , HST/ACS F814W, J, 3.6, 
4.5, 5.8, 8.0, 24 p,m. 



performed with DAOPHOT demonstrate that nearly all 
extragalactic sources are at least marginally resolved by 
IRAC, particularly at shorter wavelengths, and hence 
point source fitting is inappropriate for photometry. 

Photometry is done using the high spatial resolution 
ACS data as priors for determining an individual, object- 
specific aperture shape for extracting the IRAC fluxes. 
The advantages of this are deeper photometry and more 
accurate color measurements between low and high reso- 
lution observations. Because the IRAC beam is so large 
(compared to ACS), using the higher resolution data 
for shape information will allow us to deblend other- 
wise overlapping sources and therefore do photometry 
to fainter limits (below the confusion limit) than with 
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TABLE 1 
Observational Data 



Telescope 


Waveband 


CentralA 


Avg. 


Exptime a 


Area 


3<t Depth 




Avg. FWHM b 


Instrument 




(microns) 


(hours/pixel) 


(□') 


(AB mag.) 


(MJy) 


(arcseconds) 


i^nanara Ado-1 




z-8 rvcV 




28.0 


290 


o ci \y in — 16 c 
o.b X 10 




5.0 


r aiomar riaic Lr o 


u 


U.oD 




6.0 






U.Uo 


1 n 
l.U 


Palomar Hale LFC 


cr' 


0.47 




6.2 


450 


27.1 


0.05 


1.3 


Palomar Hale LFC 


r' 


0.63 




4.1 


450 


26.6 


0.09 


0.9 


Palomar Hale LFC 


i' 


0.77 




1.7 


450 


26.0 


0.15 


1.2 


HST ACS 


F814W 


0.81 




1.4 


260 


28.6 


0.01 


0.1 


MMT Megacam 


z' 


0.91 




3.6 


580 


25.9 


.16 


1.0 


NOAO Mayall Flamingos 


.1 


1.20 




0.5 


220 


21.4 


10.2 


1.2 


MMT SWIRC 


.1 


1.2 




5.2 


25 


24.4 


0.65 


0.8 


Palomar Hale WIRC 


J 


1.25 




0.8 


120 


21.8 


7.1 


1.1 


MMT SWIRC 


H 


1.6 




0.2 


225 


23.0 


2.4 


0.8 


Palomar Hale WIRC 


H 


1.63 




0.7 


120 


22.1 


5.4 


1.1 


Palomar Hale WIRC 


K 


2.15 




1.6 


210 


21.2 


12.3 


1.1 


Spitzer IRAC 




3.6 




8.0 


380 


25.7 


0.2 d 


1.9 


Spitzer IRAC 




1.5 




10.0 


380 


25.9 


0.17 d 


1.9 


Akari IRC 




4.5 




1.8 


100 


23.5 


1.5 


4.4 


Spitzer IRAC 




5.8 




8.0 


380 


26.3 


0.11 d 


1.9 


Spitzer IRAC 




8.0 




10.0 


380 


26.3 


0.11 d 


2.1 


Akari IRC 




11 




2.1 


100 


21.0 


15 


6.2 


Akari IRC 




15 




1.0 


100 


20.1 


33 


6.5 


Akari IRC 




18 




0.8 


100 


19.8 


15 


6.7 


Spitzer MIPS 




24 




1.0 


225 


20.8 


17.3 


5.8 


Spitzer MIPS 




70 




1.6 


160 


14.8 


3200 


18.6 



a Total and maximum exposure times are presented in [|2] 

b Average FWHM as measured from a number of point sources in the final image. 
c erg/s/cm 2 

d IRAC depths are given as measured 95% completeness limits, see 812.1.11 



IRAC alone. The disadvantage to using this technique 
is that we must assume the galaxy morphology is the 
same across the wavelength range of the high to low res- 
olution data. This assumes that ACS F814 light is co- 
spatial with IR light. At high redshifts this assumption 
will be reasonable because in those galaxies we cannot 
resolve any age or SF gradients across the galaxy. For 
low redshift galaxies we rely on the Rayleigh-Jeans side 
of the SED being less-sensitive to age so that F814 and 
IRAC are measuring the same population of stars. This 
assumption is not valid for all galaxies and must be con- 
sidered when using the catalog for science goals. The use 
of high resolution data as priors for low resolution pho- 
tomet ry has been used on similar datasets in the litera- 
ture (iFernandez-Soto et alj 119991: iPapovich et alj 120011: 
Shaplev et all 120051: iGrazian et all 120061 : lLaidler et all 
2007). 

This method begins by first running source detection 
and photometric extraction on the coadded IRAC images 
using a matched filter algorithm with image backgrounds 
determi ned using the mesh back ground estimator in SEx- 
tractor (|Bertin fc Arnou ts 1996) . This catalog is merged 
with the HST j ACS catalog. For every object in that 
catalog, if the object is detected in ACS then we use the 
ACS shape parameters to determine the elliptical aper- 
ture size for the IRAC images. ACS shape parameters 
are determined by SExtractor on isophotal object profiles 
after deblending, such that each ACS pixel can only be 
assigned to one object (or the background). For objects 
which are not detected in ACS, but which are detected 
in IRAC, we simply use the original IRAC SExtractor 
photometry. Because of the larger IRAC beam, for all 
objects in the catalog we impose a minimum semi-major 
axis of 2". In all cases aperture corrections are computed 



individually from PSFs provided by the SSC based on the 
aperture sizes and shapes used for photometry 7 . 

Final aperture photometry was performed using cus- 
tom extraction software written in IDL and based on the 
APER and MASK_ELLIPSE routines with the shape in- 
formation from SExtractor, from either ACS or IRAC 
as described above, using local backgrounds. Fluxes are 
measured in all IRAC bands at the location of all ACS 
detections without applying a flux limit (IRAC detection 
limits are discussed below). Because we use local back- 
grounds, the measured fluxes of objects near the confu- 
sion limit should have a larger scatter than those non- 
confused objects, but will on average be the correct flux. 
This will not effect the photometric redshifts, as it will 
likely shift all IRAC points up or down, but not relative 
to each other. Figure [J] shows the difference between us- 
ing the IRAC photometry alone and using the ACS based 
photometry as des cribed here. This pl ot is a direct ana- 
log of Figure 2 in lLaidler et al.l (|2007l ) which shows the 
comparison of fluxes for the software package TFIT. We 
do not see the systematic offset in fluxes that they see at 
low flux levels. 

Determining the detection limits of the IRAC data is 
complicated by varying exposure times across the field, 
source confusion, and our use of ACS locations as pri- 
ors for photometry. Field coverage is a complex function 
of spacecraft dynamics, see Figure [5] and [31 Within the 
central 180 square arcminutes of the field, the median ex- 
posure time per pixel is eight hours at 3.6 & 5.8 microns 
and ten hours at 4.5 & 8.0 microns. The maximum ex- 
posure time in the center is 28 hours and 21 hours at 3.6 
& 5.8 and 4.5 & 8.0 microns respectively. A total of 40 

7 http:/ /ssc. spitzer. caltech.edu/irac/psf.html 
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Fig. 3. — Top: The full 3.6 fxra data image. Bottom: Grey-scale 
coverage map of the 3.6 fira data with grey-scale ranging from 1 to 
1000 images, each of 100s exposure time. North is up, East is to 
the right. The actual NEP is 3.5 degrees away. 

hours of observations went into covering the entire field. 
The confusion limit in the IRAC data arises because both 
the PSF (2"FWHM) and the number of sources is large. 
We have used ACS priors to be able to deblend confused 
sources, and therefore our detection limit is below the 
nominal confusion limit. One effect of using ACS for 
deblending is that the depth of our photometry will be 
color dependent. We have IRAC photometry below the 
confusion limit for all objects which are detected in the 
F814W band, but cannot deblend objects for which their 
is no F814W counterpart. Therefore, our photometry 
is deeper for sources which have bluer spectral shapes, 
whereas it is confusion limited for the redder, IRAC only 
detections. Because of these three complexities, there is 
no one single value for the detection limit of the survey. 

We are able to measure nominal 95% completeness 
limits in the IRAC passbands from a number count dia- 
gram by measuring where the number of observed objects 
drops below 95% of the expected number of objects. The 
expected numbers are calculated by fitting a straight line 
to the brighter flux number counts and assuming that the 
faint number counts would also follow that line in the ab- 




F(IRAC) in microJy 

Fig. 4.— Difference between ACS based (F(ACS)) and IRAC- 
alone based (F(IRAC)) photometry for the 18,300 sources in com- 
mon in the catalog with flux measurements in greater than five 
wavebands. Points and error bars show the biweight center and 
dispersion of the distrubution within each bin. This allows for 
robust outlier rejection. Note that bin size changes with flux. 




Exposure Time (hours) 

Fig. 5. — Cumulative distribution of the number of pixels versus 
exposure time in hours for the central 15' diameter region. The 
solid line shows the coverage for the 3.6 & 5.8 fim images, the 
dashed line is for 4.5 & 8.0 fim images. 

sence of confusion. Figure [6] shows the number counts in 
the four IRAC bands with a best fit line to the detections 
to show where the number of measured fluxes drop below 
95% of the expected numbers. The figure also shows the 
larger confusion in the shorter wavelength bands. The 
measured and expected number counts are lower in the 
longer wavelength bands. Final values at 3.6, 4.5, 5.8, 
and 8.0 fim are 0.2, 0.17, 0.11, 0.11 /xJy respectively (see 
Table[T]). A more complete analysis of the number counts 
will be presented in a future paper. 

2.1.2. MIPS lOfim 

The 70/mi observations were taken with the Multi- 
band Imaging Photometer for Spitzer (MIPS; Riek e et all 
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Fig. 6. — Number counts in the IRAC bands showing the com- 
pleteness limits for each band. The histogram is the measured 
number of sources and the solid (red) line is a linear fit to the ris- 
ing portion of the histogram. Completeness is determined where 
the measured number of sources drops below 95% of the expected 
number of sources based on the fits. 



2004) with a total exposure time per pixel of 1.6 hours 
covering 160 square arcminutes to a limiting 4c depth 
of 4.3 mJy (Program ID 30222). The data were taken 
as a 5x2 grid of pointings (5.3' x 2.6' field of view) 
in small scale photometry mode with 10 second ex- 
posures. The raw 70/xm data were downloaded from 
the SSC archive and were reduced using the Germa- 
nium Reprocessing Tools (GeRT, version 20060415). The 
pipeline was optimized using the same parameters de- 
rived for the ultra-d eep GOODS- North photometry data 
(|Fraver et al. 2006b). We used a combination of the col- 
umn median filter and a high-pass median time filter per 
detector to remove the instrumental artifacts ("streak- 
ing"). The filtering was done in two passes to remove 
the filtering artifac ts ("negative side-lo bes" ) around the 
brightest sources (jFraver et alj l2006a|) . A calibration 
factor of 702 MJysr" 1 per MIPS-70-unit (|Gordon et all 
[2001 was assumed, and the d ata have been co l or cor - 
rected by dividing by 0.918 (|Stansberrv et alj |2007[) . 
This color-correction factor is for a constant v S v SED 
across the filter bandpass, but is appropriate to within 
2% for a wide range of galaxy SEDs (S v oc v~ a , a = 0- 
3). 

The data were combined using the SSC mosaicing soft- 
ware (mopex). PSF fitting has been found to be more 
accurate for MIPS 7 0/im data than aperture photometry 
(| Gordon et al.ll2007j ). Source extraction was done using 
the Astronomical Point-Source Extraction (apex) tools 
within MOPEX. The 70 /xm peaks were fitted using the 
Point-source Response Function (PRF; FWHM = 18'.' 6) 
imag e derived from the extragalactic First Look Sur- 
vey (|Fraver et alj I2006ah because there are not enough 
isolated bright sources in the field to accurately derive 
our own PRF. Sixty-one sources were extracted with a 



signal-to-noise ratio (SNR) greater than four (4.3mJy). 
To avoid spurious sources at 70 jtim only those 42 sources 
with 24 fim. counterparts were included in this study. We 
don't expect any sources to be detected at 70/xm that are 
not also detected at 24/iin which we confirm by eye. The 
remaining unmatched sources are blends that we are un- 
able to separate or noise peaks. We applied an aperture 
correc tion of 1.15 for e mission outside of the adopted 
PRF (jFraver et al.ll2009h . 

2.1.3. MIPS 24fim 

The 24^m data were taken with MIPS on Spitzer with 
a total exposure time of one hour per pixel covering 225 
square arcminutes to a limiting 3<7 depth of 17.3 // Jy 
(Program ID 1868). The data were taken in large-field 
photometry mode with a 30-second exposure time. A 
3x3 MIPS field of view grid was mapped and repeated 
five times, with multiple dithers and chops. There were a 
total of 1080 separate exposures with a final total depth 
of 60 minutes per pointing on the sky. The MIPS data 
were processed by the Spitzer Science Center into indi- 
vidual image BCDs (pipeline version S12.4). However, 
substantial "jailbar" artifacts, as well as a significant gra- 
dient, were visible. All of the individual pointings were 
forced to a common background by applying an additive 
constant to the entire frame. A "delta-dark" was then 
generated from the median of all frames; the great de- 
gree of dithering in the data allows this process to reject 
all actual celestial objects in the frames from the me- 
dian stack. That stack was then adjusted to a median 
overall zero value, and then subtracted from all the data. 
It currently is not known whether the gradient effect is 
additive or multiplicative, although our experience with 
other Si:As arrays of this kind strongly suggests (from a 
physical basis) that it is additive. However, we reduced 
the data both ways, and found no difference at any de- 
tectable level. The data were then coadded using the 
MOPEX software package onto exactly the same projec- 
tion system as used for IRAC, albeit with 1'.'2 pixels. 

IRAF daofind was used for object detection. We 
supply the code with the PSF FWHM (5'.'8) and back- 
ground sigma values taken by examining the image. 
Daofind then counts the flux within an annulus of di- 
ameter FWHM and flags any set of pixels as a detec- 
tion where that flux is above a threshold of five sigma. 
To deal with confused sources, we perform object detec- 
tion iteratively. After the first run through daofind, all 
objects are subtracted from the image using a PSF de- 
termined from the detected objects. Daofind is then 
re-run on the residual image. Due to the large PSF and 
the lack of nearby galaxies in the field, with the exception 
of only a handful of galaxies, all MIPS detections appear 
as point sources. Photometry on all detected sources is 
done with the IRAF task ALLSTAR which fits PSF's to 
groups of objects simultaneously. An aperture correc- 
tion of 1.4 is applied for flux beyond the 6.5 pixel radius 
at which the PSF star was normalized. This correction 
factor is calculated from a curve of growth based on the 
composite PSF star. Using this method the 3cr detection 
limit is 17.3/xJy. These noise properties are compara- 
ble to the GOODS slightly longer exposure (77 minute) 
dataset that has a 3ct limit of 12/iJy. 

2.2. Akari Space Telescope 
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The 4, 11, 15, and 18/nm data were t aken with the 
Infrared Camera (IRC; | Onaka et al"1l2007T) on the Akari 
Space Telescope (|Murakami et al.ll2007ft with a total ex- 
posure time of 1 - 2.5 hours per band per pixel covering 
100 square arcminutes to a 3a depth of 1.5, 15, 33, and 
45 /xJy respectively. The details of this data set will be 
provided in Egami et al. (2009, in preparation), so here 
we only describe some basic data parameters. 

The IRC was used to obtain images at 11, 15, and 
18 iim during 2006, October. These three Akari/IRC 
bands bridge the Spitzer imaging wavelength gap be- 
tween IRAC 8 /im and MIPS 24 ^m. A 4 fj,m image was 
also obtained since the 4 and 11 /im data were taken si- 
multaneously for the same region by design. Each image 
covers a field of view of ~10'xl0' (i.e., four times larger 
than that of IRAC) with a pixel scale of 2'/3-2'/5. The 
field center coordinates are 17:40:00 +69:00:00 (J2000). 

Low-resolution infrared spectra were also taken for 
sources in the same 10' x 10' area using the slit-less multi- 
objec t spectroscopy (MOS) mode of IRC (|Ohvama et al.1 
12007ft . When combined, obtained spectra cover a wave- 
length range of 1.8-26.5 /jm 8 with a resolving power of 
20-50. The integration times were several hours. 

2.3. Palomar and KPNO near-IR 

The J, H, a nd K s observations were taken with the 
WIRC camera ([Wilson et al.ll2003ft on the Palomar 200- 
inch Hale telescope and the FLAMINGOS camera on the 
KPNO 4m telescope. Total exposure times at Palomar 
were 0.8, 0.7, and 1.6 hours per pixel covering 120, 120, 
and 210 square arcminutes to a 3a limiting depth of 21.8, 
22.1, and 21.2 mags (AB) in J, H, and K s respectively. 
KPNO was used solely for the J band with a total ex- 
posure time of 0.5 hours per pixel covering 220 square 
arcminutes to a 3a limiting depth of 21.4 mag (AB). 
Palomar data were taken on 2006 May 5, 2008 May 15- 
18, and 2008 June 27-28 and the KPNO data were taken 
on 2006 August 29-30. In both cases the large format 
cameras are well-matched to the field size allowing effi- 
cient observing. 

The Palomar and KPNO data were reduced using a 
custom IDL pipeline (available from the authors). Blank 
sky images are generated from the actual survey data 
to serve as a temporally windowed measure of the dark 
current and sky. An initial pass through the data gen- 
erates a first-order image reduction on which SExtrac- 
tor is used to detect and mask objects and bad pixels 
within the data. A second pass utilizes the generated 
masks to reject pixels from the stacking process that cre- 
ates the blank sky images. These "sky darks" were then 
subtracted from the data. Additional tasks repair other 
image defects, such as quadrant-level floating bias un- 
certainties. The data were flat-fielded using dome flats. 
Each frame was individually calibrated both photomet- 
rically and astrometrically to 2MASS, which compen- 
sates for drifts in calibration due to non-photometric 
conditions. Individual images were re-projected using 
SWARP 9 onto a common projection center used for all 
the other datasets. The data were then coadded using 
IRAF/IMCOMBINE, with relative weighting of the im- 

8 There is a gap in coverage at 13.4-17.5 /im because the grism 
covering this range has been damaged. 

9 http:/ /terapix.iap.fr/rubrique.php?id_rubrique=49 



ages based on their noise properties. The photometric 
and astrometric calibration of the final coadded mosaic 
was then further refined by a final rematch to the 2MASS 
survey. An object catalog was extracted using SExtrac- 
tor, with a detection threshold of 2a. The SExtractor 
MAG-AUTO photometric magnitude was used. Seeing 
was typically 1.1-1.5". In general, the WIRC data are 
used for the survey as that dataset is deeper than the 
FLAMINGOS dataset. WIRC 3a detection limits vary 
from center (22.5,23.1,22.5) to the edge (21.8,22.1,21.2) 
in J, H, K (AB) respectively. The FLAMINGOS 3a de- 
tection limit is 21.4 mag (AB) in the J-band. 

2.4. MMT near-IR 

The J- and H-band observations were ta ken with the 
SAO Widefield InfraRed Camera (SWIRC: iBrown et all 
12008ft at the 6.5m MMT with a total exposure time of 7.1 
and 2.1 hours covering 25 and 225 square arcminutes to a 
3a limiting depth of 24.4 and 23.0 respectively. The data 
were taken over the course of two observing runs. The 
first set of observations were taken during the four nights 
of 2005 April 19-22. Conditions were not photometric, 
and the field was observed through thin cirrus. Total in- 
tegration times were 14 minutes at both J- and H-band 
at each of nine positions arranged in an overlapping 3x3 
configuration of a 5.1'x5.1'field of view to cover roughly 
the central 15' x 15' of the IRAC dark field. The total 
integration time was built up out of short, dithered 60s 
(J) and 30 s (H) exposures. The seeing was variable but 
remained below 1" FWHM for most of the run. The sec- 
ond set of observations were carried out during the three 
nights of 2006 May 12-14. During this campaign only the 
central 5' x 5' of the IRAC cal field (the subregion with 
the deepest overall IRAC integration) , was observed, and 
only at J-band. A total of 257 70s dithered exposures 
were acquired through variable cirrus with sub-arcsecond 
seeing for a total combined exposure time of 5 hours. 

The individual exposures were processed using stan- 
dard methods. All frames were dark- and bias-subtracted 
using dark frames taken at the beginnings and ends of the 
corresponding nights. A custom mask was constructed 
using a median stack of all the 60 s J-band exposures 
from the 2005 run, so as to exclude large regions of bad 
and noisy pixels from subsequent analysis. The celestial 
coordinates, which were initially only accurate to within 
20", were refined and verified with the wcstools pack- 
age by reference to bright 2MASS stars. The images 
were flattened with skyflats. This was an iterative pro- 
cess whereby a crude, initial mosaic was first constructed 
to reveal faint objects not easily detected in the individ- 
ual exposures. Objects in this initial mosaic were then 
identified and masked, and the deeper object masks were 
subsequently translated to the coordinate frames of all 
individual frames. In this way, with the much deeper 
object masks, much-improved skyflats were created from 
the stack of individual frames. Specifically, this tech- 
nique eliminated the negative artifacts that arose from 
repeated placement of faint-but-significant objects on the 
same parts of the detector array by the dither pattern. 

Because none of the data were taken on photometric 
nights, the photometric calibrations were referenced to 
an ensemble of well-detected 2MASS stars. Final mo- 
saics were constructed using 0'.'15 pixels by coadding the 
individual frames using sigma clipping to eliminate cos- 
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mic rays. Two wide, shallow 15' x 15' mosaics were made 
having 3a depths of 22.7 (J) and 23.0 (if) AB magni- 
tudes, based on the noise measured within a 3" diameter 
aperture. In the smaller but deeper mosaic from 2006 a 
3cr depth of 24.4 AB mag was estimated within the same 
aperture. The FWHM in the final mosaics was roughly 
0'.'75. 

Photometry was measured in the final mosaics with 
SExtractor, using MAG_AUTO as a measure of total 
magnitude, thus not requiring an aperture correction. 
For the wide-field mosaics, SExtractor was used in two- 
image mode. In this case, the IRAC 3.6 /im mosaic was 
used as the detection image, and the J- and fZ-band mo- 
saics were used as the measurement images. For the 
deep J-band mosaic, SExtractor was used in single-image 
mode. Given the 0'.'15 pixels and the $.'75 seeing, six con- 
nected pixels each individually detected at the 1.5 a level 
were required for a source detection. The resulting deep 
J-band catalog contains over 13,000 sources. 

2.5. MMT optical 

The z' observations were taken with the Megacam cam- 
era at the 6.5m MMT with a total exposure time of 3.6 
hours per pixel covering 580 square arcminutes to a 3a 
limiting depth of 25.9 mag (AB). The observations con- 
sisted of a series of 500 s exposures taken on the nights 
of 2007 June 1-4 under photometric conditions, except 
for the night of 2007 June 1, when some thin cirrus was 
present. Seeing varied from as low as ft' 6 to as high as 
1'.'7, and averaged about l'.'O. The data were reduced 
interactively using standard techniques on a dedicated 
pipeline machine hosted at SAO. Bias-subtracted and 
flattened exposures were treated to remove cosmic rays 
and bad pixels before calculating the coordinates using 
stars from the USNOB1.0 catalog and correcting the pho- 
tometry for off-axis scattered light. The resulting expo- 
sures were spatially registered to a common coordinate 
system. All frames taken on the night of 2007 June 1 were 
then flux-calibrated using exposures from the photomet- 
ric nights. All the frames were then scaled appropriately 
and coadded to create the final z' mosaic. The mosaic 
covers roughly 24' x 24' covering the entire IRAC survey 
field where the center has the deepest exposure time of 
3.6 hours tapering towards the field edge. 

The z' mosaic was photometered with SExtractor in 
single-image mode, using MAG-AUTO as a total mag- 
nitude. Because the seeing was measured to be approx- 
imately l'.'O in the combined image, the detection algo- 
rithm was configured to identify all objects having at 
least 12 connected 0'.'32 pixels lying la above the back- 
ground. Using the coverage map, portions of the mosaic 
observed three or fewer times were masked, so that pho- 
tometry was measured only for sources observed at least 
four times by Megacam. We estimate that the depth 
reached is 25.9 AB mag, 3a, within a 2"aperture. 

2.6. EST / AGS Optical 

HST I ACS F814W observations were taken with a 
total exposure time of 1.4 hours per pixel covering 260 
square arcminutes to a 3a limiting depth of 28.6 mags 
(AB) (Program ID 10521). ACS WFC is a mosaic of two 
2048 x 4096 SITE CCDs with a pixel scale of 0.05"/pixel 
(15/x pixels), so that the full field of view per exposure 



is 11.3 square arcminutes. Observations were made over 
50 orbits between November 28, 2006 and December 07, 
2006 using SIDE-2 electronics. The observations were 
divided into 25 pointings, and each pointing was in turn 
divided into eight roughly 10 minute exposures for a to- 
tal of 200 exposures. We used a small dither between 
exposures at each pointing for cosmic ray rejection and 
to cover the gap between the two ACS CCDs (2.5"). We 
also employ a ~ 0.5' overlap between pointings. The en- 
tire ACS mosaic is contained within the Spitzer IRAC 
map. 

The ACS pipeline calacs was used for basic reduction 
of the images including overscan and bias subtraction, 
cosmic ray rejection, dark subtraction, and flat field- 
ing. Beyond this base reduction we developed our own 
method of further reducing and mosaicing the 25 point- 
ings. The three largest challenges in doing this were the 
large amount of memory required to combine 200 160Mb 
fits files, getting the astrometry correct given distortions 
and relative shifts between images caused by thermal ef- 
fects in the telescope, and the variation in background 
level due to the four amplifiers used to readout the two 
CCDs, multidrizzle alone is unable to account for 
these effects. Our pipeline uses a combination of the 
IRAF tasks multidrizzle, tweakshifts, and ccmap 
along with automatic background subtraction by quad- 
rant and SWarp for the final mosaicing. The final com- 
bined mosaic image is 1.7GB in size and covers a ~ 20' 
diameter rough circle in the center of the IRAC dark field 
with the native ft'05 per pixel resolution and O'.'l seeing. 

The photometry zeropoint was taken from the head- 
ers. SExtractor was used with the same basic setup as 
the Palomar optical data for source detection and ex- 
traction (see pT7| . The SExtractor MAG_AUTO is used 
in the catalog which uses a Kron total magnitude when 
there are no near neighbors to bias the measurement and 
a corrected isophotal magnitude otherwise. The 3a de- 
tection limit for point sources is I = 28.6(^4_B). There 
are ~ 51,000 sources in the catalog with detections in 
both ACS and IRAC. 

2.7. Palomar Optical 

The u', g', r', and i' data were taken with the Large 
Format Camera at the 200" Hale telescope at Palomar 
with a total exposure time of 6.0, 6.2, 4.1, and 4.7 hours 
per pixel covering 450 square arcminutes to a 3a limit- 
ing depth of 27.2, 27.1, 26.6, and 26.0 mags (AB) re- 
spectively. Observing runs occurred on 2004 July 22-24 
and 2004 August 08-10. All observing nights had less 
than 50% moon illumination. The LFC is a mosaic of 
six 2048 x 4096 CCDs with a pixel scale of 0.18"/pixel 
(15/z pixels), so that the full field of view per exposure 
is roughly a 24' diameter circle. Images were randomly 
dithered between exposures to account for the gap be- 
tween chips (15"). Integration times were typically 360 
seconds in r' and f, 900 seconds in u', and 600 seconds 
in g'. u'-band science and calibration data were binned 
2x2 pixels. 

Reduction of this Palomar data follows standard proce- 
dures. The images were cross-talk and overscan corrected 
using IRAF tasks. A combined bias frame was then 
made from roughly five overscan corrected individual bias 
frames and applied to the data to remove the large scale 
bias. No dark current correction was required. Pixel-to- 
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pixel sensitivity variations were corrected in all images 
using 15 median-combined dome flats taken throughout 
the run with 7000-9000 counts each. Although twilight 
flats were taken during the August run, they were not 
required in the reduction, i'-band images were fringe cor- 
rected using a median smoothed, object removed, bad- 
pixel masked combination of all science frames. 

The IRAF tasks mscgetcatalog, msccmatch and 
MSCCSETWCS were used to find and apply x and y shifts 
and rotations to each frame by tying it to the USNO- 
A2 positions. Fits were done manually by ruthlessly re- 
jecting outliers while maintaining even coverage over the 
field. A combined bad pixel mask was made for all cosmic 
rays(CR), hot pixels, and satellite trails for each image 
using CRAVERAGE and masking by hand. The bad pixel 
mask was used with the IRAF task fixpix which does a 
linear interpolation over the effected pixels. To remove 
gradients in the background we use MSCSKYSUB to fit a 
second order Legendre to the background of each image 
and subtract it. SWARP was then used to re-project and 
combine together all images in each filter. 

Photometric calibration was performed in the usual 
manner using a field of SDSS well measured stars at a 
range of airmasses. Photometric nights were analyzed 
together; solutions were found in each filter for a com- 
mon magnitude zero-point. All nights except August 
10 were photometric. Those exposures taken in non- 
photometric conditions were individually tied to the pho- 
tometric data using roughly 10 stars well distributed 
around each frame. 

We use SExtractor to both find all objects in the com- 
bined frames, and to determine their shape parameters. 
The detection threshold was defined such that objects 
have a minimum of six contiguous pixels, each of which 
are greater than 1.5a above the background sky level. 
We choose these parameters as a compromise between 
detecting faint objects in high signal-to-noise regions and 
rejecting noise fluctuations in low signal-to-noise regions. 
This corresponds to 3cr point source detection limits of 
27.2, 27.1, 26.6, and 26.0 mags (AB) in u', g', r', and i' re- 
spectively. Shape parameters are determined in SExtrac- 
tor using only those pixels above the detection threshold. 

When making the catalog, we use an aperture magni- 
tude for objects fainter than 22.5 and a Kron magnitude 
from SExtractor for all objects brighter than 22.5. We 
choose a kron factor of 2.5 and a minimum radius of 3.0 
pixels as can be s hown to give the closes t estimate to 
total magnitudes (jBertin fc Arnoutslll996h . The fainter 
objects are all point sources and so do not benefit from 
shape information in their photometry. Aperture sizes 
are 3'/6, 5'.'6, 4'.'1, and 4'.'5 for u', g', r', and i' respec- 
tively. Aperture sizes were chosen based on a curve of 
growth of a bright star in the field. We choose large 
apertures (~ 4 x FWHM) to enclose 'all' of the flux and 
therefore not require an aperture correction. 

2.8. Palomar Optical Spectroscopy 

Spectroscopy was done with the COSMIC spectro- 
graph at the Palomar Hale 200" telescope. COSMIC, 
at prime focus, has a 13.6' field of view, and 0.4" pixels. 
Observations were made over four nights in June of 2007 
& 2008 with the 300 1/mm grating with a dispersion of 
3 A per pixel. We chose a slit- width of 1.5" to match 



our 1 - 1.5" seeing. The optical band covered by this 
instrument, roughly 3500 to 9000A, includes such spec- 
tral features as CaH&K, [Oil], [OIII], Ha, H/3, Ed, G 
band, and the 4000 Abreak. During both runs we were 
able to observe a total of 11 slitmasks of ~ 25 galaxies 
each with exposure times of on average 80 minutes di- 
vided into multiple exposures. One HgAr lamp and one 
flat exposure was taken through each mask at the begin- 
ning of the night for calibration. Targets were selected to 
be galaxies (non-stellar PSF's) brighter than r=21(AB) 
with priority in mask design given to those with MIPS 
24 or 70/im detections to boost the chance of seeing an 
emission line and thereby getting a secure redshift. 

Reduction was done with IRAF mainly through the 
BOGUS2006 10 scripts. Prior to running BOGUS, im- 
ages were overscan and bias subtracted. BOGUS itself 
does a 2D reduction including flat-fielding, cosmic ray 
removal, sky subtraction, fringe suppression and combi- 
nation of frames. The same reduction is performed on 
both science images and arcs. The standard IRAF tasks 
of APALL, identify, and DISPCOR were used to wave- 
length correct, trace, and extract the spectra with a sec- 
ondary background subtraction for minor level changes. 
One dimensional spectra were extracted for a total of 200 
galaxies with measurable continuum. Although spectro- 
scopic standards were observed during each run, we did 
not flux calibrate the spectra as we are only interested 
in determining redshifts (see 2]). 

2.9. Chandra A CIS -I 

Chandra ACIS-I observations were taken with a total 
exposure time of 28 hours per pixel covering 290 square 
arcminutes to a limiting depth of 3.6 x 10 _16 ergs/s/cm 2 . 
The camera is a 2 x 2 mosaic of four 1024 x 1024 frontside- 
illuminated CCDs arranged to be tangent to the mirror 
assembly in the focal plane. The pixel scale is 0.49"/pixel 
(24/x pixels), so that the full field of view per exposure is 
roughly 290 square arcminutes. Observing was done in 
three separate observations at different pointing angles 
to reduce the effect of the gap between chips (11"); 50ks 
on 2007 July 23, 22ks on 2007 September 28, and 28ks 
on 2007 September 29, for a total exposure time of lOOks 
(observation ID's 7359, 8471, 9595). Observations were 
done in timed exposure mode (TE) with very faint for- 
mat (VFAINT) in order to measure events in 5 x 5 pixel 
islands instead of smaller regions which has the advan- 
tage of better determination of background events such 
as cosmic rays. We know from the ROSAT survey that 
there are no bright objects in the field such that pile-up 
of flux would be a problem in these larger event regions. 
The data were pre-processed by the Chandra X-ray cen- 
ter with standard data processing (SDP) version 7.6.11.1. 
Further reduction was done with the Chandra Interactive 
Analysis of Observations software (CIAO 4.0). This ver- 
sion of the SDP includes identifying and correcting for 
hot pixels and cosmic ray afterglows ( residual charge on 
frames after a cosmic ray hit). 

In creating an event list we include background clean- 
ing with the code ACIS_PROCESS_events which flags pos- 
sible background events based on counts in the outer ring 
of the 5x5 event island so that they can be removed 

10 https://zwolfkinder.jpl.nasa.gov/ stern/homcpagc/ 
bogus.html 
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in the filtering phase. We then filter the event list for 
the number of pixels above a threshold and flags for bad 
pixels or cosmic rays. Secondly, frames which are not 
considered good science frames due either to the tele- 
scope motion, angle, or temperature are removed from 
the event list. Lastly we filter the data on energy range, 
keeping only energies between 0.3 and 7 KeV where the 
particle background is fairly flat and quantum efficiency 
is the highest. 

Because the observations were taken in three different 
exposures we correct the astrometry before combining 
individual observations into an exposure map and fluxed 
image. To do this we use a preliminary run of wavde- 
tect to make a source list and then REPROJECT_ASPECT 
to reproject the images to a reference image (in this case 
the first, longest observation ID 7359). MERGE_ALL is 
used to create the final combined image and exposure 
map. After combination we examine the data for back- 
ground flares and remove anomalous high background 
times from the image. We perform source detection 
and photometry on the full 0.3 - 7 Kev band data us- 
ing WAVDETECT which operates by correlating the data 
with "mexican hat" wavelet functions on size scales from 
one to eight pixels. The significance threshold parameter 
in WAVDETECT is set such that we expect approximately 
one false detection per 1000 x 1000 pixel area, or ap- 
proximately four in our field. There are 121 total source 
detections. 

Photometry with WAVDETECT sums the total counts 
within the source cell defined by the shape which encom- 
passes all positive count pixels around the source detec- 
tion. Backgrounds are determined in the wavelet analy- 
sis. To convert from count rate to flux u nits we assume 
an average spectral shape of T = 1.4 ( Giaccon i et alJ 



12002. iFd^TetaTI l200l lLaird et al.ll2009f ). Along with 
the local galactic NH value of 4 x 10 20 cm~ 2 this gives us 
energy conversion factors of 5.49 x 10~ 12 , 1.8 x 10 -11 , 
and 1.08 x 10~ n in the s oft, hard, and total b ands re- 
spectively. Using T = 1.7 ([Manners et al.ll2003h changes 
these conversion factors by ~ 10%. The minimum num- 
ber of counts we have detected in the full band is 3.3±2.0, 
corresponding to 3.6 ± 2.2 x 10~ 16 erg/cm 2 /s. 

3. THE CATALOG 

The catalog is built by matching sources from each 
band to the remainder of the catalog within a radius that 
is similar to that band's PSF size. The disadvantages of 
this method are miss-matches due to noise detections and 
deblending imperfections. A noise detection specifically 
in the ACS band will propagate through to the IRAC 
bands because of our photometry methods (see Section 
I2.1.1jl . Saturated stars with diffraction spikes often pro- 
duce spurious signals. Fortunately this field has few/no 
bright stars by design (see Section [T]). The brightest star 
in the field is 12.2 mag (AB) is F814W. The main source 
of error in the matching will be incorrect matches be- 
cause of multiple sources within the match radius. For 
this reason we match catalogs with similar PSF sizes. For 
those bands with especially large PSFs we are careful to 
match sources with other large PSF bands, and not di- 
rectly to our smallest PSF catalog. Specifically the MIPS 
24/im dataset is matched to the IRAC source list, as the 
catalog nearest in PSF size to itself. Secondly the MIPS 
70/im dataset is matched to the MIPS 24/im data. This 



method means that we will lose any MIPS 24/im sources 
which are not detected at IRAC wavelengths and simi- 
larly any 70/mi sources which are not detected at 24/im. 
We expect these objects to be extremely rare, if they 
exist at all, and therefore they will be searched for in a 
different manner. 

For this catalog we do not use matched apertures, but 
instead we calculate a "total" flux in each band. We do 
not want to lose information by smoothing all the bands 
to the same resolution, but instead keep all information 
possible at each wavelength. To this end, we use aperture 
corrections where necessary, as discussed above. 

By far, the HST ACS F814W catalog has the most 
number of sources and so dominates the spatial distri- 
bution of sources. As mentioned in 32.1.11 the ACS cat- 
alog is used as the detection image for the IRAC data. 
Other optical and near-IR bands are matched to that 
merged catalog of ACS and IRAC. In total there are 
86,815 sources in the matched catalog. This catalog will 
be made publicly available at a future date. At this time 
all of the space-based data are available from their re- 
spective archives. 

As a check of the catalog photometry, we present color- 
color diagrams of stars in the field. Stars are a good 
population on which to test our photometry since they 
follow well defined tracks in color space. We select stars 
from our ACS catalog. A plot of the SExtractor ISOAREA 
parameter as a function of magnitude clearly separates 
a ridge of unresolved point sources from the resolved 
sources because unresolved sources have smaller ISOAREA 
than galaxies at any given magnitude. Although we have 
used our highest resolution data to pick out the stars, 
there are undoubtedly interlopers from higher rcdshift 
compact galaxies or quasars. In total we have ~ 150 
point sources with detections in u', g', r', i', J, and H. 
In Figure [7] we compare the optical and near-IR colors of 
these stars to thos e from synthetic stellar sp ectra from 
the BaSeL library (|Leieune fc Sc hacrcr 20W). The com- 
parison sample includes synthetic stars with masses from 
0.4-0.8 to 120-150 M© and metallicities Z=0.0005 to 0.1 
measured at ages between 1000 yr and 16-20 Gyr. We 
find good agreement between the real star colors and 
those from the theoretical models. 

4. PHOTOMETRIC REDSHIFTS 

The combined IRAC and ACS catalog contains over 
50,000 objects which makes acquisition of spectroscopic 
rcdshifts impractical. Even confirmation spectroscopy of 
red galaxies at z — 1 in our three candidate clusters 
will require many nights on 8-10m class telescopes and 
is therefore also impractical. In lieu of spectroscopy we 
use our extensive multi-wavelength, broad-band catalog 
to build spectral energy distributions (SEDs) using up to 
13 bands (u\ g', r', i\ F814W, z', J, H, K, 3.6, 4.5, 5.8, 
8.0/im) from which we derive photometric redshifts. 

These SEDs are fit with template spectra from 
IPolletta et alJ (|2007l ) 11 . These templates have been used 
successfully by a num ber of surveys at a range of redshifts 
for a l l galaxy types (lAdami et al.l 120081: INegrello et all 
I2009t iSalvato et all l2009t lllbert et all 12009ft . Because 
the templates include the IR spectral region (empiri- 
cally built from observations in SWIRE), we find them 

11 http:/ /www. iasf-milano.inaf.it/~polletta/templates/swire_templatcs. html 
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Fig. 7. — Optical and near-IR color color diagrams as a 
check of photometr y . Col ors from synthetic stellar spectra from 
Lcjeune & Schacrcr (2001) for a set of stars covering a range in 
mass and metallicity are shown as black points. The 150 stars 
from the dark field are shown as triangles. 

the best choice to use as models for this dataset. We 
use 15 templates including ellipticals, spirals, star form- 
ing galaxies, and AGN. Photometric redshifts are calcu- 
lated using Hyperz; a chi-squared minimization fitting 
program including a co rrection for interstellar reddening 
( Bol zonella et all 120001 : iCalzetti et all 12000). An exam- 
ple of a well fit galaxy is shown in Figure [5] Errors in 
photometric redshifts are determined by comparing the 
photometric redshifts with spectroscopic redshifts. 

Spectroscopic redshifts were determined using both 
IRAF tasks emsao and xcsao. Specifically emsao 
searches the spectrum for both absorption and emission 
lines which it correlates with a given line list. XCSAO 
cross-correlates the spectrum with known galaxy tem- 
plates which allows us to use features like the 4000 A 
break and the rest of the spectral shape to identify red- 
shifts. Both techniques were used together to arrive at 
the best fit redshift for each galaxy. We used 17 spec- 
tral templates of galaxies and AGN from the compilation 




Fig. 8. — Example spiral galaxy with a well fit photometric red- 
shift by hyperz. The thumbnails show a 12" X 12" box in r' , HST / 
ACS F814W, K, 3.6, 4.5, 5.8, 8.0, 24 fim, and the fitted SED. The 
title of the SED plot gives the best fit photometric redshift, the 
probability that that redshift is good, and the name of the best-fit 
template. Triangles are data used in the fit, the line is the fit, and 
the cross is the 24 (im data point, not used in the hyperz fit, but 
plotted here for completeness. 

of th e HST Calibration D atabase System (iFrancis et al.l 
[19911: iKinnev et ail H991 ; ICaTzetti et al.lll994D . We do 
not use the SWIRE SED templates here (as we do for the 
photometric redshifts above) because the optical spec- 
tra in the SWIRE templates is based on modeling and 
their true strength lies in their IR treatment. Since our 
spectroscopy is in the optical, we use the Calibration 
Database System templates because they are based more 
reliably on UV and optical observations. We applied a 
very strict requirement that all emission and absorption 
features in the ID spectra were confirmed by eye on the 
2D spectra and that multiple lines be identified in all 
cases to avoid incorrect redshift determination due to 
cosmic rays or noise features from sky line subtraction. 

We were able to successfully determine redshifts for 
87 galaxies. This represents a conservative sample of 
'good' redshift determinations defined to have either high 
signal-to-noise emission lines or multiple absorption lines 
and good cross correlations. We then compare the spec- 
troscopic to photometric redshifts to quantify the error 
on the latter (Figure [9]) . There are cases where Hyperz 
has failed to fit the correct redshift. In these galaxies 
it appears that Hyperz has not chosen a template that 
goes through any of the photometric points, completely 
missing the SED shape. This could be due to a fail- 
ure of hyperz or the errors in the photometry. Those 
galaxies, as characterized by a y 2 value greater than 50, 
are not included in this comparison sample. The error 
on the photometric redshifts is 0.064(1 + z). Note that 
this error is quoted as a function of redshift and so takes 
into account the increasing scatter with z. This accuracy 
is similar to other IRAC based multi- wavelength s tudies 
(jBrodwin et al.ll200a iRowan-Robinson et al.ll2008f) . 

5. INITIAL SCIENCE RESULTS 

Given this unique dataset, we are carrying out sci- 
entific programs to exploit the available information. 
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Fig. 9. — Comparison of spcctroscopically and photometrically 
determined redshifts for 87 galaxies. The scatter implies an error 
on the photometric redshifts of 0.064(1 + z). 

iKrick et "ail (|2008f ) have discovered a large scale structure 
at redshift one consisting of three galaxy clusters. These 
clusters show evidence of a build-up of the red cluster se- 
quence at redshift one. In a study of the star formation 
and AGN fraction in those clusters, IKrick et al.l (|2009D 
have found evidence for star formation in star forming 
galaxies to look the same regardless of environment, but 
a dense environment does strengthen the suppression of 
star formation and AGN activity in clusters over the field 
from redshift one to the present. 

Ongoing work continues on finding and character- 
izing high redshift galaxies [z > 6), compton thick 
AGN, the lowest m ass brown dwarfs, population III SNe 
(|Frost et al.l I2009H - the reddest objects in the universe 
from their 24/im — r colors, and other interesting phe- 
nomenon, the study of which is uniquely suited to this 
field. As this is the only mid-IR dataset with time do- 
main information on a 5+ year time-line, we are explor- 
ing the variabl e and transient objects in the field, from 
stars to AGN dHund et alj|2006t ISurace et al.ll2007l and 
Foesneau et al. 2009, in preparation). One such appli- 
cation is to study the extragalactic variability in mid-IR 
to set constraints on the nature of dust around AGN. 

As a representation of the data, we include here a dis- 
cussion of IRAC color color diagrams to bot h differenti- 
ate between AGN and s tar forming galaxies (jLacv et al.l 
120041: iStern et all [2005) and demonstrate the effective- 
ness of our ACS based IRAC photometry. Fig \W\ shows 
all sources originally detected in the four IRAC bands 
before adding ACS shape information. Here we have 
separated the sample based on SED classifications into a 
spiral, elliptical, AGN, and stellar sample. Those SEDs 
best fit with composite templates are not colored in the 
plot. Deep surveys like the IRAC dark field, in op- 
position to the shallower, wide-area, surveys like the 
Syitzer first look survey and the IRAC shallow survey 
(|Eisenhardt et al.l l2004h show a larger amount of con- 
tamination in the AGN wedges fro m intermediate red- 
shift spiral and ell i ptical galaxies (iSaiina" eTaTI [20051 
iBarmbv et all [200l iBarger et all l2008f) . Those AGN 
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Fig. 10. — IR AC color color diagram in the fashion of a) 
ILacv et al.l J2004T) . b) IStern et all f?0(m . Sources arc color coded 
by best fit SED shape. Green, yellow, bl ue, and red i ndicat e stars, 
ellipticals, spirals, and AGN. While the ILacv et al.l (120041) wedge 
has high completeness for the AGN sam ple, there is significant 
contamination, and the Stern et al. (2005) wedge shows lower com- 
pleteness but less contamination. 

which are blue- ward (left) of the lLacv et alJ (|2004f ) wedge 
are AGN at redshift ~ 1.0 that have a stellar peak in ad- 
dition to a power law shape, and that stellar peak has 
been shifted into the IRAC 3.6 /zm band. We see from 
these plots that our SED fitting produces nice results for 
the morphologies of the galaxies. 

Secondly we discuss the effect of our ACS based IRAC 
photometry on the scatter in the color colo r diagrams. 
We d o see larger scatte r than the original ILacv et ail 
(|2004h and IStern et all (|2005l) diagrams because this 
is a deeper survey. The larger scatter is also seen in 
other deep suryeys li ke that in the Extended Groth Strip 
(jBarmbv et al.ll200l . We use this scatter in the color 
color diagrams to quantify the difference between using 
the IRAC only photometry and the ACS based IRAC 
photometry. We use the same sample of ~ 11, 000 galax- 
ies originally detected in all four IRAC bands to make 
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color-color distributions of both IRAC only and ACS 
based IRAC photometry. The change in scatter is quan- 
tified with the F-test where the ratio of variances gives 
probabilities that the scatters in the distributions are the 
same. Those probabilities are less than 1 x 10~ 14 indicat- 
ing that the scatter is statistically significantly reduced if 
we use the ACS based photometry instead of the IRAC 
only photometry. 
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